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Abstract 


Nomenclature 

basis  vector 

direction  cosine  array 

or  damping  matrix 

unit  vector 

or  eliminated 

current  frame  of  reference 

geometric 

general  constraint  relationship  function 

inertial  frame  of  reference 

stiffness  matrix 

linear  operator 

mass  matrix 

number  of 
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The  initial  version  of  the  General  Rotorcraft  Aerome- 
chamcal  Stability  Program  (GRASP)  was  developed  for 
analysis  of  rotorcraft  in  steady,  axial  flight  and  ground 
contact  conditions.  In  these  flight  regimes,  the  mate¬ 
rial  continua  of  the  rotorcraft  may  experience  deformations 
which  arc  independent  of  time.  GRASP  can  obtain  this 
steady-state  solution  and  can  solve  the  eigenproblem  associ¬ 
ated  with  perturbations  about  such  a  steady-state  solution. 
GRASP  it  the  first  program  implementing  a  new  method^ 
ogy  for  dynamic  analysis  of  structures,  parts  of  which  may 
be  experiencing  discrete  motion  relative  to  other  parts.  Ap¬ 
plication  of  this  new  methodolagfrto  GRASP,  including  sub¬ 
structuring,  frames  of  reference,  nodes,  finite  elements  and 
constraints,  is  described  in  the  paper.  GRASP  combines  fea¬ 
tures  usually  found  in  traditional' finite  element  programs 
with  those  usually  found  in  the  multibody  programs  used 
for  spacecraft  analysis.  The  structure  is  decomposed  into  a 
hierarchy  of  substructures,  and  discrete  relative  motion  be¬ 
tween  substructures  is  treated  exactly.  Deformation  of  con¬ 
tinua  is  treated  by  the  finite  element  method.  The  library  of 
finite  elements  includes  a  powerful  nonlinear  beam  element 
that  incorporates  aeroelastic  effects  based  on  a  simple  non¬ 
linear,  aerodynamic  theory  with  unsteady  induced  inflow.*. 
The  analytical  bases  for  the  aeroelastic  beam  element  and  ] 
the  screw  constraint  are  presented,  and  the  important  role ' 
of  geometric  stiffness  in  the  formulation  is  illustrated. 


generalised  force  ^  : 

generalised  coordinate 
position  * 

retained  I 

superordinate  (parent)  frame  of  reference 
the  identity  matrix 
variation  of 
virtual  work 
virtual  rotation 
Levi-Civita  symbol 
magnitude  of  Euler  rotation 
or  components  of  a  small  rotation 
Euler- Rodrigues  parameters 


Convent  ions 


=  the  Gibbsian  vector  named  a 
=  the  unit  vector  a 
=  the  steady-state  value  of  a 
—  the  perturbation  value  of  a 
s  the  associated  skew  symmetric  matrix  for  the 
vector  a 

=  a  after  steady-state  deformation 
*  a  after  steady-state  deformation  and 
perturbation 

=  the  a  associated  with  &  relative  to  c 
=  the  transpose  of  the  matrix  a 
=  the  ith  component  of  the  vector  a  in  the  b  basis 
=  the  ith  component  of  the  vector  a 
or  the  ith  instance  of  the  quantity  a 
=  vector  cross  product  of  a  and  b 


Introduction 


Over  the  years,  research  analytes  that  treat  various  as¬ 
pects  of  rotorcraft  aeroelastic  stability  have  evolved  from 
very  simple  models1  that  were  used  to  gain  insight  into 
the  physical  phenomena,  into  more  complex  models*-4  en¬ 
deavoring  to  accurately  represent  realistic  rotorcraft  config¬ 
urations.  At  the  same  time,  large  simulation  programs?-* 
which  include  many  sophisticated  features  not  present  in 
the  research  codes,  were  being  developed  to  solve  a  variety 
of  rotorcraft  problems.  These  programs  made  significant 
contributions  to  the  understanding  of  rotorcraft  modeling. 
As  discussed  in  Ref.  9,  however,  they  were  not  designed 
with  rotorcraft  stability  as  their  primary  application,  nor 
do  they  possess  the  depth,  generality,  and  modeling  flexi¬ 
bility  needed  to  handle  the  full  range  of  rotorcraft  stability 
problems. 
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It  is  against  this  background  that  the  General  Rotor- 
craft  Aeromechanical  Stability  Program  (GRASP)  was  de¬ 
veloped.  The  initial  version  of  GRASP  was  designed  specif- 
ieally;fae  the  purpose  of  calculating  the  aeromechanical  sta¬ 
bility  ot  rotorcraft  in  axial  flight  and  ground  contact.  Con- 
•  figurations  or  conditions  that  would  result  in  periodic  or 
otherwise  time-dependent  motion  cannot  currently  be  ac¬ 
commodated.  GRASP  is  designed  to  provide  a  library  of 
solution  methods;  the  initial  version  provides  a  steady-state 
solution  dn$  an  eigensolution.  The  steady-state  solution  is 
<5ne  in  wHKh  the  deformations  within  each  material  con- 
tjhuum  of  the  structure  remain  constant,  although  there 
■lay  be  discrete  motions  between  the  material  continua. 

eigtnproblem  solution  is  one  in  which  the  deformation 
‘  in  each  material  continuum  is  a  constant,  reference  de¬ 
lation  (typically  the  steady-state  solution)  plus  a  per¬ 
turbation.  Eigenvalues  and  eigenvectors  can  be  obtained 
from  the  linearized,  perturbation  equations.  Expansion  of 
GRASP  capabilities  to  other  flight  regimes  can  be  accom¬ 
plished  by  adding  solutions  to  the  library.  A  description  of 
the  solution  procedures  lies  outside  the  scope  of  this  paper. 
However,  a  detailed  discussion  may  be  found  in  Ref.  10. 

GRASP  is  the  first  program  which  has  been  written 
utilizing  the  new  methodology  described  in  a  companion 
paper?1  Research  in  spacecraft  dynamics  provided  the  ba¬ 
sis  for  this  methodology,  which  incorporates  body  flexibility 
with  the  large  discrete  motions  previously  available  with 
multibody  programs?1  The  concept  underlying  the  method¬ 
ology  is  that  a  structure  is  modeled'as  composed  of  substruc¬ 
tures  which  may  undergo  discrete  motions  relative  to  one 
another.  The  substructures  themselves  are  deformable  elas¬ 
tic  continua  which  are  represented  by  finite  element  models 
in  GRASP.  Thus,  through  the  use  of  libraries  of  nodes,  ele¬ 
ments,  and  constraints,  GRASP  provides  the  powerful  mod¬ 
eling  capabilities  normally  available  with  the  finite  element 
method.  In  addition  GRASP  has  the  capability  for  allowing 
large  discrete  motions  between  substructures. 

The  primary  intent  of  this  paper  "is  to  describe  the  im¬ 
plementation  of  the  new  methodology  in  GRASP.  In  so  do¬ 
ing,  some  of  the  more  important  features  of  the  theoretical 
basis  for  GRASP  are  described.  A  complete  presentation 
of  the  equations  for  all  aspects  of  the  GRASP  analysis  is 
beyond  the  scope  of  this  paper  but  may  be  found  in  the 
GRASP  Theoretics^  Manual10 

The  first  section  of  this  paper  describes  the  GRASP 
implementation  of  the  decomposition  of  the  structure  into 
substructures.  The  second  and  third  sections,  respectively, 
explain  the  concepts  of  substructure  reference  frames  and 
nodes.  The  fourth  section  describes  the  implementation  of 
finite  elements  in  GRASP.  Three  finite  elements  are  cur¬ 
rently  in  the  GRASP  library:  the  aeroelastic  beam,  the  air 
mass,  and  the  rigid-body  mass.  A  description  of  the  aero¬ 
elastic  beam  is  presented  as  an  example  of  the  way  a  finite 
element  is  implemented  in  GRASP.  The  fifth  section  outlines 
the  development  of  constraints.  At  an  example  of  how  the 
constraint  equations  are  derived,  a  kinematics!  development 
for  the  screw  constraint  is  presented.  Through  this  exam¬ 
ple,  the  role  that  geometric  stiffness  plays  in  the  constraint 
equations  is  also  demonstrated. 


Subatruct  uring 

The  first  step  of  the  methodology  is  the  decomposition 
of  the  structure  into  substructures  such  that  discrete  motion 
occurs  only  between  substructures  (e.y.,  the  rotor  and  fuse¬ 
lage  of  a  rotorcraft  must  be  in  separate  substructures).  In 
GRASP,  substructuring  is  extended  to  a  hierarchical  collec¬ 
tion  of  substructures.  That  is,  the  structure  is  decomposed 
into  substructures,  which  in  turn  are  decomposed  into  (sub- 
Substructures  (a  sub-substructure  will  just  be  called  a  sub¬ 
structure)  and  so  on.  The  process  is  continued  until  all 
the  lowest  level  substructures  consist  of  a  single  finite  ele¬ 
ment.  The  user’s  input  to  GRASP  uses  a  slightly  different 
terminology.  The  elements  are  referred  to  as  element-type 
subsystems  and  all  of  the  other  substructures,  including  the 
complete  structure,  are  referred  to  as  system-type  subsys¬ 
tems.  Discrete  motion  is  prohibited  within  element-type 
subsystems  but  may  be  specified  within  a  system-type  sub¬ 
system. 

There  is  a  great  deal  of  flexibility  in  the  way  a  structure 
can  be  decomposed  into  a  hierarchy  of  substructures,  pro¬ 
viding  a  powerful  organisational  tool  for  modeling.  Besides 
isolating  discrete  motions,  some  of  the  uses  for  substructur¬ 
ing  are  to  separate  physical  components,  to  provide  logical 
separation  by  function,  to  isolate  changing  components  for 
parametric  studies,  and  to  allow  natural  coordinatisation  of 
components. 

The  hierarchical  organization  of  the  substructures  is 
represented  in  GRASP  as  a  tree  information  structure.  The 
root  of  the  tree  is  created  internally  by  GRASP  and  corre¬ 
sponds  to  the  degrees  of  freedom  of  the  complete  structure 
after  applying  constraints.  The  only  child  of  the  root  is  the 
complete  structure  (a  system-type  subsystem).  From  that 
point,  the  tree  branches  out  into  the  substructures  which 
branch  out  into  substructures  until,  at  the  most  detailed 
level,  called  the  leaves  of  the  tree,  each  substructure  con¬ 
tains  a  single  finite  element  (element-type  subsystem).  The 
GRASP  user  conveys  the  position  of  each  substructure  in 
the  tree  by  including  the  name  of  its  parent  subsystem  in 
the  definition  of  the  subsystem  containing  the  substructure. 
It  should  be  emphasized  that  a  substructure  is  a  subset  of 
its  parent,  not  a  branch  off  of  it.  Thus,  the  topology  of  the 
structure  is  not  restricted  to  a  tree  configuration. 

Within  GRASP,  instead  of  associating  the  term  sub¬ 
structure  with  its  physical  embodiment,  the  term  may  be 
used  to  identify  a  collection  of  abstractions  or  data  struc¬ 
tures  associated  with  the  model  of  the  substructure.  From 
that  point  of  view,  a  substructure  is  thought  of  as  being 
composed  of  a  frame  of  reference,  a  set  of  nodes,  a  set  of  sub¬ 
structures,  a  set  of  finite  elements  (which  are  viewed  inter-  • 
nally  as  just  more  substructures),  and  the  set  of  constraints 
which  can  be  used  to  extract  the  independent  generalised 
coordinates  for  the  substructure  from  all  of  the  generalised 
coordinates  for  the  substructure. 

Each  substructure  is  also  associated  with  a  state  vector,  - 
a  residual  vector,  and  a  set  of  mass,  damping,  and  stiffness 
matrices.  Starting  at  the  lowest  level  (a  substructure  con¬ 
sisting  of  a  single  finite  element),  the  elements  of  these  arrays  , 
are  associated  with  the  substructure’s  generalised  coordi¬ 
nates.  This  includes  the  generalized  coordinates  for  the  ele¬ 
ment’s  nodes,  for  the  element’s  nodeless  variables,  and,  since 1 


an  dement  is  a  substructure  and  thus  has  a  frame  of  refer¬ 
ence,  for  the  element’s  frame  generalised  coordinates.  Af¬ 
ter  applying  the  constraints  for  the  substructure  (element), 
only  the  independent  generalised  coordinates  remain.  The 
independent  generalised  coordinates  for  the  substructure  be¬ 
come  part  of  the  complete  set  of  generalised  coordinates  in 
the  parent  structure.  This  process  continues  with  each  par¬ 
ent  structure  containing  the  generalized  coordinates  of  its 
children.  Finally,  at  the  root,  the  elements  of  the  arrays  are 
associated  with  every  independent  generalized  coordinate  in 
the  whole  structure  (after  applying  the  set  of  constraints  as¬ 
sociated  with  the  whole  structure). 

The  general  concept  of  hierarchical  substructuring  is 
that  everything  that  belongs  to  a  substructure  is  just  a  sub¬ 
set  of  what  belongs  to  the  parent.  However,  in  GRASP,  even 
though  the  node  in  a  child  is  logically  part  of  the  parent, 
only  the  independent  generalized  coordinates  from  a  node 
are  available  in  the  parent.  The  node  itself  is  not  available 
in  the  parent.  This  poses  a  difficulty  when,  for  instance, 
constraining  a  node  in  a  child  to  a  node  in  the  parent.  It 
is  desirable  to  have  the  constraint  defined  within  one  sub¬ 
structure  for  a  variety  of  programming  reasons  and  because 
the  constraint  is  naturally  written  in  one  frame  of  reference. 
In  GRASP,  this  is  resolved  by  creating  another  node  in  the 
parent,  associated  with  the  same  physical  region  of  the  ma¬ 
terial  continuum  as  the  node  in  the  child.  Since  there  are 
now  two  nodes  and  two  sets  of  generalized  coordinates  as¬ 
sociated  with  the  same  set  of  physical  degrees  of  freedom, 
a  constraint  is  required.  Such  a  node  promotion  constraint 
(there  are  also  node  demotion  constraints)  requires  that  the 
two  nodes  behave  identically.  In  GRASP,  these  constraints 
are  entirely  invisible  from  the  user’s  point  of  view.  GRASP 
automatically  generates  nodes  and  constraints  for  its  inter¬ 
nal  use.  The  user  merely  specifies  that  nodes  in  different 
substructures  are  to  be  constrained  together. 


Reference  Frames 

In  GRASP,  the  overall  frame  of  reference  for  the  prob¬ 
lem  is  inertial.  However  the  frame  of  reference  for  the  struc¬ 
ture  as  a  whole  may  have  a  nominal  motion  which  is  a  spec¬ 
ified  constant  velocity,  acceleration,  and  angular  velocity 
relative  to  the  inertial  frame  of  reference.  Spin-up  is  not 
currently  supported  in  GRASP.  The  frame  for  the  structure 
as  a  whole  is  taken  to  be  coincident  with  the  inertial  frame 
initially  (since  the  position  and  orientation  of  the  inertial 
frame  are  not  relevant  to  solving  the  problem). 

In  addition,  there  is  a  reference  frame  associated  with 
each  substructure  contained  in  the  model  hierarchy.  These 
substructure  frames  of  reference  are  selected  to  participate 
in  the  gross  motion  of  the  substructures.  During  calcula¬ 
tions  for  a  substructure,  the  inertial  position  and  orienta¬ 
tion  of  the  superordinate  (parent)  frame  and  the  position 
and  orientation  of  the  current  frame  relative  to  the  superor¬ 
dinate  frame  are  available.  Given  these,  the  inertial  position 
and  orientation  of  the  current  frame  can  be  determined  by 
using  the  appropriate  chain  rules  as 


Rrt  =  Rrs  +  Rsl 
CFI  =  cFscst 


(1) 


Each  frame  of  reference,  including  the  one  for  the  whole 
structure,  introduces  six  redundant  generalized  coordinates 
associated  with  the  rigid-body  motions  of  the  frame.  For  the 
steady-state  problem,  the  translational  generalized  coordi¬ 
nates  are  the  measure  numbers  for  the  steady-state  position 
of  the  frame  relative  to  the  reference  position  of  the  frame 
expressed  in  the  basis  associated  with  the  steady-state  po¬ 
sition  of  the  frame,  R^,r .  The  rotational  generalized  co¬ 
ordinates  are  the  Euler- Rodrigues  parameters  (also  called 
Rodrigues  parameters  or  components  of  the  Rodrigues  vec¬ 
tor)  for  the  steady-state  orientation  of  the  frame  relative  to 
the  reference  orientation  of  the  frame,  df  F  - 

GRASP’s  definition  of  Euler-Rodrigues  parameters  dif¬ 
fers  from  the  usual  definition11  by  a  factor  of  2,  viz., 

d?A  =  2eA.tan?  (2) 

where  a  rotation  of  magnitude  9  occurs  about  a  unit  vector 
i  (an  Euler  rotation)  which  carries  the  frame  A  into  coinci¬ 
dence  with  frame  B.  The  direction  cosine  matrix  associated 
with  an  arbitrary  set  of  those  Euler-Rodrigues  parameters 
arranged  as  a  column  matrix  denoted  by  d  >* 

rBA  (1  —  <frT<t>/4)&  +  M>T ~  4  ... 

C  “  1+^/4  (3) 

The  tilde  symbol  is  defined  by 

A<y  =  -  <iji>Aa  (4) 

(Recall  that  the  Levi-Civita  symbol,  *,>*,  is  1  if  i,j,k  is  an 
even  permutation  of  1, 2, 3;  -l  if  i,j,  k  is  an  odd  permutation 
of  1,2,3;  and  0  otherwise.)  Besides  the  simplicity  of  the 
conversion  to  and  from  direction  cosines,  Euler-Rodrigues 
parameters  have  the  advantage  of  being  free  of  singularities 
for  rotations  up  to  180? 

The  generalised  coordinates  for  the  eigenproblem  rep¬ 
resent  a  perturbation  beyond  the  steady-state  values.  The 
translational  generalized  coordinates  are  hp»r' .  Since  the 
perturbations  are  infinitesimal,  the  rotational  generalized 
coordinates  can  be  associated  with  the  measure  numbers  of 

a  rotation  vector  .  The  direction  cosine  matrix  asso- 
r  * 

ciated  with  these  generalized  coordinates  is 

Cr“r  =  A-&.  (5) 

To  eliminate  the  redundant  generalized  coordinates  and 
establish  the  position  and  motion  of  the  frame  relative  to  the 
substructure,  a  constraint  must  be  specified.  In  GRASP, 
the  constraint  library  currently  contains  three  frame  con¬ 
straints:  the  fixed  frame,  the  periodic  frame,  and  the  rotat¬ 
ing  frame  constraints.  The  prescribed  constraint,  normally 
used  to  constrain  nodal  generalised  coordinates,  can  also 
be  used  to  prescribe  specific  frame  generalised  coordinates. 
Each  of  the  frame  constraints  specifies  the  motion  of  the 
frame  in  terms  of  the  superordinate  (parent)  frame. 

When  GRASP  begins  the  calculations  associated  with 
a  substructure,  the  first  thing  that  is  done  is  to  calculate 
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the  position  sad  orientation  of  the  frame  of  reference  rela¬ 
tive  to  the  fuperordinate  frame  and  the  steady-state  iner¬ 
tial  motion  of  the  frame  of  reference.  The  relative  position 
and  orientation,  B?3'  and  CF'3‘ ,  may  be  generated  from 
the  frame  constraint  parameters  and  the  frame  generalized 
coordinates  for  the  current  and  superordinate  frames.  The 
steady-state  inertial  motion  of  the  current  frame  of  reference 
can  be  obtained  from  the  same  information  as 


in  the  steady-state  problem  only  (and  is  constrained  to  be 
zero  in  the  eigenproblem).  qA  and  qf  are  the  cyclic  velocity 
gradients  which  are  both  active  only  in  the  eigenproblem 
(and  are  constrained  to  be  zero  in  the  steady-state  prob¬ 
lem).  For  the  eigenproblem,  since  the  generalised  coordi¬ 
nates  are  velocities  and  spatial  gradients  thereof,  they  sire 
actually  considered  to  be  the  time  derivatives  of  generalized 
coordinates. 


£','  =  as'/xfl''5'+i:s'/  (6) 

^'-A^x  (o*'x 


To  obtain  the  contribution  of  frame  generalized  coordi¬ 
nates  to  the  virtual  work,  virtual  displacements  of  the  frame 
generalized  coordinates  are  required.  The  virtual  displace¬ 
ments  are  just  variations  of  the  steady-state  problem  trans¬ 
lational  and  the  eigenproblem  translational  generalized  co¬ 
ordinates.  The  virtual  work  associated  with  the  steady-state 
problem  rotational  and  the  eigenproblem  rotational  gener¬ 
alised  coordinates  is  treated  differently.  The  virtual  work  is 
the  dot  product  of  the  moment  and  a  virtual  rotation  vector, 
the  measure  numbers  of  which  are  tipF\F  and  SipF',',F'  for  the 
steady-state  problem  and  eigenproblem,  respectively.  More 
detail  may  be  found  in  Ref.  10  or  the  related  discussion  in 
a  companion  paper-1 1 


Nodes 

Substructure  equations  of  motion  in  GRASP  are  devel¬ 
oped  based  on  a  finite  element  formulation.  Rather  than 
solving  for  the  field  itself,  nodes  are  used  to  define  general¬ 
ized  coordinates  that  approximate  the  field  in  the  vicinity 
of  the  node.  In  GRASP,  there  are  two  fields,  the  structural 
deformation  field  and  the  air  flow  field.  The  library  of  nodes 
supports  both  fields,  containing  one  node  for  each  field:  the 
air  node  and  the  structural  node.  Air  nodes  to  define  the 
air  flow  fields  in  the  aeroelastic  beam  and  air  mass  elements. 
Structural  nodes  define  the  displacement  fields  in  the  aero¬ 
elastic  beam  and  rigid-body  mass  elements. 

Air  Node 

The  only  air  flow  field  currently  supported  in  GRASP 
is  the  axisymmetric  flow  field  of  a  helicopter  rotor  includ¬ 
ing  dynamic  inflow  effects.14  The  air  node  must  be  inertially 
fixed  and  located  on  the  axis  of  symmetry.  The  basis  associ¬ 
ated  with  the  air  node,  A,  has  the  4a,  axis  pointed  upwind 
along  the  axis  of  symmetry  (vertical  in  hovering  flight).  The 
inertial  velocity  of  the  air  flow,  U,  at  any  point,  P,  in  the 
flow  field  is 

npt  -  -  {it  +  V  +  +  <V)  i.,  (7) 

where  JtFA  is  the  perpendicular  distance  to  the  axis  of  the 
flow  field 

*1-  -  V<«if)’+W)’  <»> 

The  four  qA  are  the  air  node  generalized  coordinates. 
qf  is  the  axial  flow  velocity  which  is  an  active  generalised 
coordinate  for  both  the  steady-state  problem  and  the  eigen¬ 
problem.  qA  it  the  radial  velocity  gradient  which  is  active 


Structural  Node 

GRASP  is  designed  to  support  displacement  fields 
which  are  characterized  by  small  strains  but  large  displace¬ 
ments  and  rotations.  DispltLcements  are  represented  by 
the  traditional  three  translational  generalized  coordinates. 
However,  the  usual  approach  of  taking  small  rotations  about 
the  basis  vectors  as  generalized  coordinates  is  inappropriate 
for  representing  large  rotations.  Instead,  the  three  general¬ 
ised  coordinates  are  taken  to  be  the  Euler- Rodrigues  param¬ 
eters  which  specify  the  orientation  of  the  node  after  defor¬ 
mation.  Thus  the  generalized  coordinates  for  a  structural 
node  are  the  ones  that  would  be  appropriate  for  specify¬ 
ing  the  motion  of  an  infinitesimal,  massless  rigid  body  that 
occupies  the  material  point  associated  with  the  node. 

The  generalized  coordinates  for  a  structural  node  are, 
thus,  analogous  to  the  generalized  coordinates  for  a  frame 
of  reference.  The  translational  generalized  coordinates  for 
the  steady-state  problem  and  eigenproblem  are  R^tN  and 
The  rotational  generalised  coordinates  for  the 
steady-state  problem  and  the  eigenproblem  are  ^ N  and 
For  the  translational  generalized  coordinates',  the 
virtual  displacements  are  just  the  variation  of  the  transla¬ 
tional  generalised  coordinates.  For  the  rotational  general¬ 
ised  coordinates,  the  virtual  rotations  S+ and  Sip %"I'r 
are  used  in  the  steady-state  problem  and  eigenproblem. 

Elements 

In  its  current  version,  GRASP  includes  three  elements: 
the  aeroelastic  beam,  the  air  mass,  and  the  rigid-body  mass. 
The  equations  for  each  of  these  elements  were  derived  by  us¬ 
ing  the  principle  of  virtual  work.  After  a  general  description 
of  how  one  goes  about  deriving  the  equations  for  a  general 
element,  the  derivation  of  the  equations  for  a  specific  ele¬ 
ment,  the  aeroelastic  beam,  will  be  discussed. 

Element  Equations 

In  general,  the  derivation  of  the  equation  for  a  finite 
element  involves  first  selecting  a  set  of  generalized  coordi¬ 
nates  that  discretize  the  internal  displacements  of  the  ele¬ 
ment.  Since  the  element  displacements  are  also  represented 
by  nodal  coordinates,  a  transformation  between  the  nodal 
coordinates  and  element  generalised  coordinates  must  be  de¬ 
fined.  In  the  case  of  the  rigid-body  mass  element,  these  two 
steps  are  trivial,  since  the  nodal  coordinates  are  identical  to 
the  element  generalized  coordinates. 

Next,  expressions  for  the  virtual  work  of  internal,  body, 
and  surface  forces  are  derived.  The  internal  forces  arise  from 
deformations  of  the  element  itself,  while  the  body  and  sur¬ 
face  forces  are  imposed  on  the  element  form  outside  sources. 
The  rigid-body  mass  element,  for  example,  is  subject  only 
to  inertial  and  gravitational  body  forces. 
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Asroelastic  Beam  Element 

A*  an  example  of  the  derivation  of  the  equation*  for 
an  element,  consider  the  aeroelastic  beam.  It  is  the  most 
general  element  currently  included  in  the  GRASP  element 
library,  and  is  used  to  model  slender  beams  undergoing  small 
strains  and  large  rotations.  In  the  derivation,  shear  deforma¬ 
tion  and  inplane  distortion  of  the  cross  section  are  ignored. 
While  warping  rigidity  is  also  ignored,  all  other  effects  of 
warping  have  been  retained.  The  nonlinear  beam  kinemat¬ 
ics  for  this  type  of  element  have  been  formulated  and  applied 
to  the  dynamic  analysis  of  a  pretwisted,  rotating,  beam  el¬ 
ement  in  Ref.  15. 

lu  GRASP,  the  generalized  coordinates  of  the  aeroelas¬ 
tic  beam  element,  Fig.  1,  are  carried  by  a  frame  of  reference 
F,  two  structural  nodes  R  and  T,  and  sin  air  node  A.  The 
axes  of  both  the  frame  and  the  root  structural  node  are  cen¬ 
tered  at  the  shear  center  of  the  root-end  cross  section,  and 
oriented  along  the  cross-section  principal  axes.  Similarly, 
the  tip  structural-node  axe*  are  centered  at  the  shear  center 
of  the  tip-end  cross  section,  and  oriented  along  the  cross- 
section  principal  axes.  The  air  node  axes  are  centered  at . 
the  flow  field  axis  of  symmetry,  and  oriented  such  the  1-axis 
is  along  the  axis  of  symmetry. 


perturbed  for  the  dynamics,  depend  on  the  nodal  coordi¬ 
nates.  Geometric  stiffness  is  also  discussed  below  in  refer¬ 
ence  to  constraints.  Further  information  and  details  of  this 
calculation  may  be  found  in  Ref.  10. 

Elasticity 

The  elastic  equations  for  the  beam  are  based  on  the 
variation  of  the  strain  energy 


Spatial  Discretization 

The  interior  displacements  of  the  beam  are  represented 
by  four  function*  of  the  axial  coordinate  xj:  u,  and  8i. 
Bending  is  described  by  ui  and  uj,  axial  displacement  by 
«s,  and  torsion  by  0s .  These  functions  are  discretised  in 
terms  of  standard  cubic  and  linear  polynomials  so  that  the 
generalised  coordinates  at  the  root  and  tip  of  the  beam  can 
be  related  to  the  nodal  displacements  and  rotations.  In 
addition,  however,  there  are  also  generalized  coordinates, 
termed  internal  degrees  of  freedom,  associated  with  higher- 
order  polynomials. 

The  variable*  and  0j  are  expanded  in  a  set  of  polyno¬ 
mials  based  on  Ref.  16.  The  “CO"  functions  (us  and  0s)  are 
expanded  in  terms  of  0i(x)  where  x  =  x5/f.  All  functions 
after  the  first  two  standard  linear  functions  are  orthonor- 
malised.  The  “Cl"  function*  (u0)  are  expanded  in  terms  of 
/?,(*)•  Similarly,  all  functions  after  the  first  four  standard 
cubic  functions  are  orthonormalized. 

Transformation  of  Coordinates 

Since  the  nodes  have  a  different  set  of  generalized  co¬ 
ordinates  than  those  of  the  beam  element,  it  is  necessary  to 
calculate  the  beam  generalized  coordinates  in  terms  of  the 
nodal  displacements  and  rotational  variables  at  both  the 
root  and  tip  of  the  beam.  In  general,  the  beam  generalized 
coordinates  are  functions  of  the  nodal  displacements,  and 
the  orientation  of  the  tip  node  relative  to  the  root  node. 

In  addition  to  the  transformation  of  the  nodal  coor¬ 
dinates  to  beam  generalized  coordinates,  the  generalized 
forces  associated  the  beam  generalized  coordinates  must  be 
transformed  into  the  forces  and  moments  at  the  root  and 
tip  nodes.  In  the  steady-state  problem,  the  transformation 
is  straightforward,  but  tedious.  However,  in  the  dynam¬ 
ics  the  transformation  contributes  additional  geometric  stiff¬ 
ness  term*  at  both  the  root  and  tip  of  the  beam.  These  term* 
arise  because  the  transformation  equations,  which  must  also 


and  where  0(xj)  is  the  pretwist  angle,  with  0(0)  =  0,  and 
(  )'  =  d(  )/dx j.  The  components  of  the  curvature  vector,  a,, 
depend  nonlinearly  on  u'a  and  0j ,  and  the  extension  of  the 
elastic  axis,  tjj,  depends  nonlinearly  on  u'.  After  integrat¬ 
ing  over  the  cross  sectional  area,  we  obtain  the  variation  of 
the  strain  energy  as  follows: 


SU  =  J  (Fj^jj  +  M\6k\  +  M}6k2  +  (11) 

where  Fj,  Mi,  Mj,  and  M}  are  the  stress  resultants,  which 
can  be  expressed  in  terms  of  cH  and  fc<.  Details  of  this 
derivation  can  be  found  in  Ref.  15. 

Inertial/ Gravitational 

The  generalized  forces  due  to  the  motion  of  the  aero- 
elastic  beam  relative  to  an  inertial  frame  are  also  derived 
following  Ref.  15.  This  portion  of  the  derivation,  which  ig¬ 
nores  warping,  is  based  on  the  work  done  by  inertial  and 
gravitational  forces  acting  through  a  virtual  displacement. 
The  virtual  work  is  then 

SW  =  J\sup  ■EP  +  6JP-  KP)d*s 

+  6*P  ■  j‘  EPdxi  +  6±p  J\kP  +  BPP  x  £ P)d*s 

P  #P  ’  (12) 

where  £  and  Af  are,  respectively,  the  inertial  and  gravi¬ 
tational  forces  and  moments  associated  with  a  generic  point 
P  along  the  beam  axis,  RPF  is  the  position  of  that  point 
relative  to  the  element  frame,  Sup  and  6j>p  are  the  virtual 
displacement  and  virtual  rotation  at  that  point,  and  S]tF 
and  S£p  are  the  virtual  displacement  and  virtual  rotation 
of  the  frame. 


Aerodynamics 

The  aerodynamic  force*  acting  on  the  aeroelaatic  beam 
element  are  determined  from  a  quasi-steady  adaptation  of 
Greenberg'*  thin -airfoil  theory  (Ref.  17).  Two  new  *et*  of 
coordinate  axe*  muit  be  introduced  for  the  purpose  of  defin¬ 
ing  the  direction*  in  which  the  aerodynamic  force*  and  pitch¬ 
ing  moment  act.  In  Fig.  2,  the  Z  axes  are  shown  to  be  a 
set  of  dextral  axes  associated  with  the  sero-lift  line  for  the 
airfoil  section.  The  other  set  of  axes  is  the  so-called  wind 
axes  W.  For  these  axes  the  base  vector  kw,  >*  identical  to 
.  The  base  vector  is  along  the  relative  wind  vector 
(in  the  direction  of  drag)  and  is  in  the  direction  of  lift. 
As  with  other  axes,  these  axes  convect  with  the  local  beam 
cross  section  during  deformation. 

The  aerodynamic  forces  and  pitching  moment  act  at  the 
aerodynamic  center  Q  (which  is  assumed  to  be  the  quarter 
chord).  The  offset  position  of  Q  relative  to  the  origin  of 
the  local  principal  axes  P  is  Rz^kz,-  The  dynamic  wind 
velocity  vector  at  the  aerodynamic  center  denoted  by  W ^ 
is  calculated  by  subtracting  the  inertial  structural  velocity 
at  Q"  from  the  inertial  air  velocity  at  Q". 

The  angle  of  attack  is  determined  from 


w?: 

tana  =  — A- 


(13) 


The  local  air  flow  velocity  gradient  is  expressed  as 

aw?: 


gq"  = 


dR' 


(H) 


The  relative  wind  velocity  magnitude  is  given  by 


wr  =  V/W),  + W)’  <l5> 

Both  the  relative  wind  velocity  and  the  angle  of  attack  are 
time-dependent  quantities  that  depend  on  the  kinematical 
variables  for  the  aeroelaatic  beam  element,  including  frame 
motion  and  induced  inflow  velocity  degrees  of  freedom.  The 
static  and  dynamic  perturbations  of  these  quantities  are  de¬ 
termined  and  used  in  GRASP  in  their  exact  form  to  calcu¬ 
late  the  generalised  forces. 

The  applied  force  on  the  beam  is  assumed  to  be 


£  *  £*fctr"  +  +  £nt  t>Z[' 

and  the  applied  moment  is 

M  =  M'hz.. 


(16) 


(17) 


The  governing  equations  for  the  components  of  the  aero¬ 
dynamic  forces  and  moments  are  adapted  from  thin  airfoil 
theory  with  arbitrary  lift,  drag  and  pitching  moment  coeffi¬ 
cients  ci,c*,  and  em. 


Ce  =  \p.W\c,  +  ^WG%, 

V  =  \p.W\c4 

M  =  \p.W*c*cm  -  (l VG%,f  (18) 


With  the  assumption  that  the  air  exerts  a  force  on  the  struc¬ 
ture  that  is  equal  and  opposite  to  that  of  the  structure  on 
the  air,  the  virtual  work  for  a  dynamic  perturbation  can  be 
formed  as 

6W  =  l  (-SSf^Fz" +6?%lM)dxt  (19) 

where  the  6S%.<  end  6 T^„  are  the  virtual  displacements 
and  rotation  of  an  element  of  air  relative  to  the  structure, 
respectively.  These  virtual  displacements  and  rotations  in¬ 
clude  all  of  the  kinematical  variables  for  the  beam  element 
m  well  as  degrees  of  freedom  representing  frame  motion  and 
dynamic  inflow  velocity  and  velocity  gradients.  The  virtual 
work  per  unit  of  beam  element  length  done  by  the  aerody¬ 
namic  forces  and  pitching  moment  can  then  be  put  into  the 
form  of  steady-state  generalised  forces  and  linear  coefficient 
matrices  for  generalised  accelerations,  velocities,  and  dis¬ 
placements.  Details  of  calculating  the  aerodynamic  contri¬ 
bution  to  the  elements  of  M,  £7,  K,  and  Q  for  the  disefetised 
system  and  other  aspects  of  the  formulation  can  be  found 
in  Ref.  10.  It  is  worthwhile  to  note  that  the  aerodynamic 
contribution  to  the  matrix  M  is  symmetric. 

Constraints 

The  GRASP  constraint  library  contains  constraints 
serving  a  variety  of  purposes.  These  include  frame  con¬ 
straints,  element  connectivity  constraints,  nodal  and  gener¬ 
alised  coordinate  constraints,  and  constraints  generated  by 
GRASP  that  remain  internal  to  the  program  (e.y.,  the  node 
promotion/demotion  constraints  mentioned  in  the  substruc¬ 
turing  section).  All  of  the  constraints  which  are  currently  in 
the  library  have  certain  features  in  common.  This  permits 
a  general  description  of  the  way  in  which  constraints  are 
treated.  The  treatment  for  a  specific  constraint,  the  screw, 
is  then  presented  a*  an  illustration.  The  screw  constraint 
highlights  the  important  role  of  geometric  stiffness  in  dy¬ 
namical  analyses  with  nonlinear  constraints.  Further  detail 
on  the  screw  and  other  constraints  can  be  found  in  Ref.  10. 

Treatment  of  Constraints  in  General 

A  constraint  creates  dependency  among  the  generalised 
coordinates.  In  GRASP,  the  constraints  and  the  depen¬ 
dency  amongst  the  generalised  coordinates  are  eliminated 
to  that  the  actual  matrix  eigenproblem  and  the  iteration  for 
a  steady-state  solution  are  performed  only  on  the  indepen¬ 
dent  problem  degrees  of  freedom.  For  all  of  the  constraints 
currently  in  the  GRASP  library,  the  constraint  relationship 
can  be  written  in  the  special  form 


WWW? 


«*=*(«•■,.■  ••»««■„»)»  (»  =  1,...,JV*)  (20) 

Thus  the  generalised  coordinate*  related  to  the  constraint 
can  be  partitioned  into  two  sets.  And,  the  set  to  be  elim¬ 
inated  can  be  obtained  directly  from  the  constraint  func¬ 
tions  which  only  depend  on  the  set  to  be  retained.  All  of 
a  substructure's  generalised  coordinate*  can  be  obtained, 
using  the  constraint  relationships,  to  obtain  the  eliminated 
generalised  coordinate*  from  the  retained  generalised  coor¬ 
dinates.  (It  was  noted  in  the  substructuring  section  that  in 
the  case  of  the  steady-state  problem  state  vector,  the  inde¬ 
pendent  degrees  of  freedom  of  the  substructure  are  stored 
in  the  parent  substructure.) 

The  virtual  work  for  the  generalized  coordinates  asso¬ 
ciated  with  the  constraint  is 

N'  N' 

W  (21) 

i-t  1=1 

The  sum  of  the  generalized  forces  associated  with  a  gener¬ 
alised  coordinate,  Q,  may  differ  from  zero  for  two  reasons. 
First,  while  seeking  a  steady-state  solution,  equilibrium  is 
not  satisfied  and  these  residuals  ate  a  measure  of  the  error 
in  the  approximate  solution.  Second,  even  if  the  system  is 
in  equilibrium,  individual  substructures  will  not  be.  For  in¬ 
stance,  the  root  node  for  an  aeroelastic  beam  element  will 
not  be  in  equilibrium  at  the  element  level  substructure;  it 
will  only  be  in  equilibrium  in  a  higher  level  substructure 
where  the  contributions  from  all  the  elements  connected  to 
that  node  have  been  included. 

Taking  the  variation  of  Eq.  (20)  yields 

**•<  =  Y  lr^  (••  =  !,...,*•)  (22) 

;>1  "r» 

Substituting  Eq.  (22)  into  Eq.  (21)  yields 

«*-£**,  (23) 


(<  =  !,...,*•) 


Expansion  of  the  generalized  force,  Q,  including  both  elim¬ 
inated  and  retained  terms,  yields 

N 

Q.  =  Q.  +  («  =  1,-..,W)  (26) 

]=i 

where  the  linear  operator,  L>  contains  the  terms  normally 
associated  with  the  mass,  damping  and  stiffness  matrices, 
— ~  C ~  K ■  (The  minus  signs  are  present  because 
the  generalized  force  is  generally  regarded  as  positive  on 
the  right  hand  side  of  the  dynamical  equation,  whereas  the 
linear  coefficient  matrices  are  regarded  as  positive  on  the 
left  hand  side.) 

GRASP  calculates  the  M,  C,  and  K  matrices  for  a  sub¬ 
structure  by  adding  the  contributions  of  each  of  its  children. 
The  rows  and  columns  of  the  child  substructure’s  matrices 
correspond  to  all  of  the  generalized  coordinates  of  the  child. 
The  constraints  are  used  to  eliminate  dependent  generalized 
coordinates,  resulting  in  matrices  whose  rows  and  columns 
correspond  to  only  the  retained  generalized  coordinates  of 
the  child.  The  matrices  elements  are  then  added  to  the  ele¬ 
ments  of  the  parent  substructure's  matrices  that  correspond 
to  the  child’s  independent  degrees  of  freedom.  The  required 
transformations  can  be  found  using  the  virtual  work  for  the 
substructure. 

An  expression  of  virtual  work  expanded  about  the 
steady-state  value  may  be  obtained  by  substituting  Eq.  (25) 
and  the  eliminated  and  retained  subsets  of  Eq.  (26)  into  the 
virtual  work  expression  in  Eq.  (21).  After  discarding  terms 
of  second  or  higher  order,  the  result  is  a  constant  term  and 
two  first  order  terms  in  q.  The  constant,  first  term  is  the 
same  as  Eq.  (23)  evaluated  at  the  steady-state  equilibrium. 

The  linear,  second  term  is 


This  relationship  is  used  by  GRASP  to  incorporate  the  con¬ 
tributions  of  the  generalized  forces  associated  with  the  elim¬ 
inated  generalized  coordinates  into  the  retained  generalized 
forces.  During  calculation  of  steady-state  residuals,  the 
residual*  associated  with  the  eliminated  generalized  coor¬ 
dinates  are  transformed  and  added  to  the  appropriate  resid¬ 
uals  in  the  parent  substructure’s  residual  vector. 

The  treatment  of  constraints  for  the  eigenproblem  is  a 
little  more  involved.  For  the  eigenproblem,  the  generalized 
coordinates  are  assumed  to  be  the  sum  of  a  steady-state 
value  and  an  infinitesimal  perturbation  (i.e.,  q  =  q  +  q). 
Equations  (20),  (22),  and  the  generalized  forces,  Q  can  all 
be  expanded  in  a  Taylor  series  about  the  steady-state  value. 
Noting  that  Eq.  (20)  is  valid  when  q  =  q,  expansion  of 
Eq.  (20)  yields 


Expansion  of  Eq.  (22)  yields 


(«  =  1 . N ’)  (24) 


-£l>. 


(27) 


The  matrix  Ka  represents  the  geometric  stiffness  associ¬ 
ated  with  the  constraint.  During  assembly  of  the  matrices 
for  the  parent  substructure,  GRASP  calculates  this  geomet¬ 
ric  stiffness  and  adds  it  to  the  stiffness  matrix  in  the  parent 
substructure.  This  extremely  important  term  is  often  over¬ 
looked.  For  instance,  a  pendulum,  modeled  as  a  rigid-body 
mass  constrained  to  rotate  about  an  offset  axis  (using  a 
screw  constraint)  derives  a/1  of  its  stiffness  from  this  geo¬ 
metric  stiffness  term. 

The  linear,  third  term  is 


Y  Y  (Lr>k  +  Y  (28) 

jm  I  S.|  \  »1  / 
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Separating  the  turn  on  k  into  partial  sums  involving  only 
eliminated  or  retained  generalised  coordinates,  and  substi¬ 
tuting  Eq.  (24)  into  Eq.  (28)  for  the  eliminated  perturbation 
coordinates  yields  for  the  third  term 


AT  AT  /  N*  -s -  W  -3 - 

jmi  S.I  \  i>l  U*r>  1-1  ,r* 

.ff^T  "5T\  . 


The  quantity  within  the  parentheses  in  Eq.  (29)  can  be 
thought  of  as  defining  a  new  set  of  M,  C ,  and  K  matri¬ 
ces  in  terms  of  the  retained  and  eliminated  portions  of  the 
original  matrices.  GRASP  calculates  the  new  matrices  and 
adds  their  elements  to  the  elements  of  the  parent  substruc¬ 
ture’s  matrices. 

The  complete  definition  of  a  constraint,  then,  follows 
from  specification  of  g.  An  understanding  of  the  effect  of 
the  constraint  in  the  dynamics,  however,  cannot  be  obtained 
without  explicit  calculation  of  the  matrix  |j  and  the  geo¬ 
metric  stiffness  matrix  Ka. 


Example  Constraint:  The  Screw 

For  the  screw  constraint,  consider  two  nodes  connected 
by  a  mechanism  that  permits  translation  and  rotation. 
Specifically,  one  node  is  permitted  to  displace  relative  to 
the  other  along  an  axis,  which  is  denoted  as  the  screw  axis, 
that  is  fixed  in  the  deflected  and  rotated  coordinate  system 
of  both  nodes.  Also,  one  node  is  permitted  to  rotate  relative 
to  the  other  about  that  same  axis.  Denote  the  dependent 
node  (whose  generalized  coordinates  are  to  be  eliminated) 
by  D  and  the  independent  node  (whose  generalized  coordi¬ 
nates  are  to  be  retained)  by  I.  The  I  generalized  coordinates 
are  known  and  it  is  necessary  to  write  a  relationship  between 
them  and  the  D  .generalized  coordinates.  The  screw  axis  is 
identified  by  a  unit  vector  {,tf.  To  Simplify  the  algebra, 
two  intermediate  nodes  are  introduced  denoted  by  5  and 
M  (for  “stationary”  and  “moving”),  located  on  the  screw 
axis.  I  and  S  are  locked  together  during  deformation  of 
the  substructures  to  which  they  are  attached  as  are  D  and 
M.  M  and  S  initially  coincide  positionally  but  may  differ 
in  orientation.  Rf1,  CDt ,  RpM,  and  Rf1  are  assumed  to 
be  given. 

5lslici 

For  the  static  problem,  the  relationship  governing  the 
generalised  coordinates  and  the  contribution  of  the  force  and 
moment  acting  at  D*  to  those  at  I'  must  be  written.  The 
basic  displacement  and  orientation  relationships  govern  the 
position  and  orientation  of  the  dependent  node  in  terms  of 
those  for  the  independent  node  and  are 

b*d  =  a** ?r  +  an 

+  fi/s  +  (30) 

CCO  _  QVW  (jlttlit  (jiffs'  QS'r  qV  lplD 


where  M'  denotes  a  node  whose  position  and  orientation 
relative  to  S'  is  the  same  as  that  of  M  relative  to  5.  In  order 
to  obtain  the  component  form  of  the  kinematical  relation 
corresponding  to  the  function  g  above,  the  first  of  Eqs.  (30) 
are  written  in  component  form.  Let  *  =  »t“  where 
u  is  the  screw  displacement  and  RMS  =  0.  In  the  D  basis, 
the  position  becomes 

R%°  =  CDlCir  Crs" Cs  **  Ca 14  CM  D  Rp.1* 

-  R%M  +  CDlCtr Cl  * tie‘s"'  (31) 

+  CD/(C/r  Rf, r  +  Rrt 1  -  Rf1) 

The  equations  simplify  somewhat  if  CIS  =  Cvs"  = 
CD'M'  _  cDM  =  A.  is  nicely  expressed  as  an  Eu¬ 

ler  rotation  given  0,  the  screw  rotation,  and  the  unit  vector 
about  which  the  rotation  occurs  is  e'sc,r  =  eJJT  ,  the  screw 
axis  unit  vector.  The  second  of  Eqs.  (30)  is  implemented  in 
GRASP  so  that  first  each  of  the  direction  cosine  matrices  on 
the  right  hand  side  is  determined.  All  are  given  in  terms  of 
geometrical  parameters  for  the  screw  except  CM  M  which 
is  known  in  terms  of  the  screw  rotation  angle  0  and  Cl  1 
which  is  known  in  terms  of  the  Euler- Rodrigues  parameters 
for  the  node  I.  The  Euler- Rodrigues  parameters  at  D  can 
then  be  calculated  once  CD  °  is  known.  To  write  out  this 
nonlinear  transformation  in  detail  would  be  quite  difficult, 
but  part  of  the  philosophy  GRASP  is  that  such  equations 
are  not  to  be  written  out  in  explicit  form.  This  reduces  the 
chance  for  error  and  makes  the  code  more  readable. 

For  the  purpose  of  transforming  the  generalized  forces, 
which  for  this  constraint  are  components  of  the  actual  force 
and  moment  at  V  and  the  screw  force  and  moment,  the  vir¬ 
tual  displacement  and  virtual  rotation  components  must  be 
determined.  The  virtual  displacement  components  are  the 
same  as  the  Lagrangian  variation  of  the  vector  components, 
namely 


6R°'°  /  -  (Rf  4-  tie',"’  +  *fr)*V 


+  t\"'  0ti+eYr‘  R?'M'  60] 


The  virtual  rotation  components  are  given  by 


=  CD,(S*fl  +  e'r'60)  (33) 


A  force  and  moment  at  D'  will  contribute  to  the  virtual 
work  at  the  screw  connection  (through  6ti  and  60)  and  at  /’ 
(through  6Rf 1  and  St^j1).  The  virtual  work  done  by  f0” 
and  »* 


6W  =  (6Rf,)TFf/  +  ’fiR? rFf>'  +  iif) 

+  6u(e'"‘  )T  Ff*'  (34) 

+  60{er')T(R?"F?'  +  Aff) 


Clearly,  the  force  components  at  /’  are  just  those  at  D' 
transformed  into  the  /  basis.  The  moment  components  at 
however,  include  not  only  the  moment  at  D'  transformed 
into  the  I  basis,  but  also,  as  expected,  the  moment  of  the 
force  at  D'  about  /’  in  the  /  basis.  The  screw  force  is 
the  component  of  the  force  at  D'  along  the  screw  axis  for 


(40) 


the  deformed  structure.  Similarly,  the  screw  moment  is  the 
component  of  the  moment  at  D*  along  the  screw  axis,  aug¬ 
mented  by  the  moment  of  the  force  at  D'  about  the  screw 
axis. 

The  §J  matrix  turns  out  to  be 


CDt  -CDIkfr  CDIeYr‘  -CDIA?“'e  ft 
0  CDl  0  CDIefT‘ 


(35) 


where  the  columns  of  the  matrix  are  associated  with  the 
variations  of  the  generalized  coordinates  SRj  1 ,  Su, 

and  69  and  the  rows  of  the  matrix  are  associated  with  6Rp  D 
and  6i/>qD. 


Dynamic* 

For  the  dynamics  problem,  a  similar  relationship  gov¬ 
erning  the  generalized  coordinates  is  used  to  find  the  matri¬ 
ces  and  Ka.  Consider  the  nodes  and  the  screw  axis  in 
their  dynamic  states,  an  infinitesimal  perturbation  from  the 
static  position  and  orientation.  The  basic  equations  for  the 
position  are  similar  to  those  of  the  static  case 


&D"D'  =R°" 

+  R*'1"  +  &rr  +  &‘‘D' 


(36) 


For  the  orientation,  as  an  alternative  to  the  above,  the  vir¬ 
tual  rotations  can  be  used  directly. - 


6il> 


D"D‘ 


=  6rt> 


D"  to" 


+  6il> 


m"  to" 


+  6il> 


to  "s" 


+  6t/>s‘  r  +  s  +  6if> 


/"S' 


.ro‘ 


(37) 


Here  the  first,  third,  fourth,  and  sixth  terms  are  zero.  Pro¬ 
ceeding  as  above  and  noting  that 


Crr  =  A  -  C'irc"'  =  Crt(  A  -  9fl')Cir 

«Crr--Cr^A-#frwfVr 

6Cm"u"  =  ~69e‘ff"CM"!C,'‘ 

»  (38) 

the  j*  matrix  turns  out  to  be  the  same  as  Eq.  (35).  It 
should  be  noted,  however,  that  with  a  different  choice  of 
coordinate  bases  for  the  nodes,  such  as  is  used  for  the  frames 
in  GRASP}*  the  matrix  |j  may  not  be  the  same  for  statics 
as  it  is  in  dynamics.  The  geometric  stiffness  matrix  Ka 
turns  out  to  be 


0 

0 

0 

0 


0 

-Ff}'Rf'r 

-(«r’)T*f 

-<•}•''  )»■*/>' 


o  o 

Ppt'r  -Ff  Rf'^'eY*' 

0  0 

n  >rf}"' 

U  «•"'  Ff>‘ 

(39) 


where  the  columns  correspond  to  the  perturbations  of  the 
generalised  coordinates  r ,  6/  r ,  u,  end  8  and  the  rows 
correspond  to  the  variations  of  the  generalized  coordinates 
6Rf  1 ,  6i/rf  f ,  6u,  and  6&. 


Note  that  the  geometric  stiffness  matrix  is  not  nec¬ 
essarily  symmetric.  In  particular,  note  the  three- by-three 
rotation-rotation  term  associated  with  the  independent 
node.  It  is 


rvr 

The  reasons  for  this  are  the  choices  of  generalized  coordi¬ 
nates  (infinitesimal  rotations)  and  virtual  generalized  co¬ 
ordinates  (virtual  rotations)  for  the  equations.  A  similar 
result  was  obtained  by  Roberson  and  Likins}*  It  is  easily 
shown  that,  for  example,  use  of  Euler- Rodrigues  parame¬ 
ters  and  their  variations  to  express  both  quantities,  respec¬ 
tively,  leads  to  a  symmetric  geometric  stiffness  matrix.  The 
present  method  was  implemented  in  GRASP  because  rows 
of  the  matrix  equations  corresponding  to  the  virtual  rota¬ 
tions  are  physical  moments  whereas  the  generalized  forces 
associated  with  the  variation  of  Euler- Rodrigues  parameters 
have  no  easily  identifiable  physical  significance. 

Just  from  Eqs.  (39)  and  (40)  it  is  evident  that  geometric 
stiffness  is  an  important  phenomenon.  It  must  be  accounted 
for  under  all  circumstances  in  which  there  is  a  coordinate 
transformation  that  in  nonlinear.  Most  frequently,  this  phe¬ 
nomenon  is  associated  with  finite  rotation  of  the  coordinate 
axes  that  are  involved  in  the  transformation. 

Concluding  Remarks 

In  response  to  the  limitations  of  previous  methods 
for  analyzing  rotorcraft  stability,  GRASP  has  been  devel¬ 
oped.  The  first  version  of  GRASP  provides  both  a  nonlinear 
steady-state  solution  and  linearized  eigensolution  for  rotor- 
craft  in  axial  flight  and  ground  contact  conditions.  Time- 
varying  solutions  are  not  provided  in  the  current  version. 
GRASP  is  the  first  program  to  be  written  utilizing  the 
new  methodology  described  in  a  companion  paper? r  The 
methodology  permits  both  discrete  motions  between  por¬ 
tions  of  the  structure  and  deformation  of  the  structure. 
It  thus  combines  the  flexibility  in  structural  representation 
available  with  the  traditional  finite  element  approach  with 
the  ability  to  handle  large  discrete  motions  between  portions 
of  the  structure  available  with  the  multibody  approach  usu¬ 
ally  associated  with  spacecraft  applications. 

In  GRASP,  the  first  step  of  the  methodology,  the  de¬ 
composition  of  the  structure  into  substructures,  has  received 
additional  emphasis.  A  hierarchical  decomposition  of  the 
structure  is  supported  and  continues  to  the  finite  element 
level. 

The  second  step  of  the  methodology,  associating  a 
frame  of  reference  with  each  substructure,  is  supported  by 
a  library  of  frame  constraints.  There  must  be  a  frame  con¬ 
straint  to  specify  the  behavior  of  the  6  redundant  general¬ 
ized  coordinates  introduced  by  the  frame.  The  library  cur¬ 
rently  supports  frames  at  a  constant  relative  displacement 
and  orientation,  frames  rotating  at  a  constant  angular  ve¬ 
locity,  and  frames  replicated  periodically  about  an  axis.  By 
adding  frame  constraints  to  the  library,  additional  ways  of 
associating  a  frame  of  reference  with  a  substructure  can  be 
incorporated  in  GRASP. 

The  third  step  of  the  methodology,  developing  the  equa¬ 
tions  of  motion,  is  based  on  the  finite  element  method  in 
GRASP.  GRASP  contains  libraries  of  nodes,  finite  elements, 
and  constraints  specifying  the  connectivity  of  the  elements. 
Currently,  there  are  two  types  of  nodes:  structural  and  air; 
and  three  types  of  elements  and  associated  connectivity  con¬ 
straints:  aeroelastic  beam,  air  mass,  and  rigid-body  mass. 


By  adding  nodes  and  element*  to  the  libraries,  additional 
modeling  capability  can  be  incorporated  in  GRASP. 

Two  characteristics  of  the  GRASP  finite  element  mod¬ 
els  are  noteworthy.  First,  they  include  terms  reflecting  the 
motion  of  the  frame  of  reference  for  the  substructure  con¬ 
taining  the  element.  Second  they  include  nonlinear  effects 
and  support  deformations  which  result  in  large  displace¬ 
ments  and  rotations.  The  aeroelastic  beam  element  was 
used  to  provide  an  example  of  the  implementation  of  a  fi¬ 
nite  element  in  GRASP.  It  provides  a  powerful  modeling 
tool,  allowing  the  user  to  select  the  order  of  the  polynomial 
approximation  of  the  displacement  field  and  incorporating 
aeroelastic  effects. 

The  fourth  and  final  step  in  the  methodology,  the  spec¬ 
ification  of  constraints,  is  supported  by  a  library  of  con¬ 
straints.  In  addition  to  the  previously  mentioned  frame 
and  connectivity  constraints,  there  are  a  collection  of  con¬ 
straints  which  constrain  nodal  generalized  coordinates.  The 
library  currently  supports  prescribing  a  specified  general¬ 
ised  coordinate,  rigidly  connecting  two  nodes,  and  allowing 
screw  motion  (translation  along  and  rotation  about  a  single 
axis).  By  adding  constraints  to  the  library,  additional  ways 
of  constraining  the  motion  of  nodes  can  be  incorporated  in 
GRASP. 

One  of  the  most  powerful  constraints  is  the  screw  con¬ 
straint,  which  can  be  used  to  build  up  a  variety  of  gimbal 
or  hinged  connections.  Besides  providing  an  example  of  the 
implementation  of  a  constraint  in  GRASP,  it  demonstrates 
the  importance  of  nonlinear  effects  resulting  in  geometric 
stiffness  in  the  model. 

The  algorithms  used  for  solving  the  equations  obtained 
with  the  methodology  were  outlined  in  a  companion  paper?1 
Solutions  obtained  with  the  methodology  have  been  excel¬ 
lent.  Comparisons  of  GRASP  results  with  experimental  re¬ 
sults  in  the  nonlinear  analysis  of  a  cantilever  beam  are  re¬ 
ported  on  in  Refs.  9  and  15  and  in  a  companion  paper?* 
Ref.  20  reports  on  the  correlation  of  GRASP  results  with 
experimental  data  for  a  torsionally-soft  rotor. 
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